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SAMPLING OF REAL MULTIVARIATE POLYNOMIALS AND 
PLURIPOTENTIAL THEORY 


ROBERT J. BERMAN AND JOAQUIM ORTEGA-CERDA 


Abstract. We consider the problem of stable sampling of multivariate real 
polynomials of large degree in a general framework where the polynomials are 
defined on an affine real algebraic variety M, equipped with a weighted mea¬ 
sure. In particular, this framework contains the well-known setting of trigono¬ 
metric polynomials (when M is a torus equipped with its invariant measure), 
where the limit of large degree corresponds to a high frequency limit, as well as 
the classical setting of one-variable orthogonal algebraic polynomials (when M 
is the real line equipped with a suitable measure), where the sampling nodes 
can be seen as generalizations of the zeros of the corresponding orthogonal 
polynomials. It is shown that a necessary condition for sampling, in the gen¬ 
eral setting, is that the asymptotic density of the sampling points is greater 
than the density of the corresponding weighted equilibrium measure of M, as 
defined in pluripotential theory. This result thus generalizes the well-known 
Landau type results for sampling on the torus, where the corresponding critical 
density corresponds to the Nyqvist rate, as well as the classical result saying 
that the zeros of orthogonal polynomials become equidistributed with respect 
to the logarithmic equilibrium measure, as the degree tends to infinity. 


1. Introduction 

1.1. Background. By the classical Whittaker-Shannon-Kotelnikov sampling the¬ 
orem a band-limited signal / on the real line M, normalized so that its frequency 
is in [—1,1] may be recovered from its values at the points tj = jir where j ranges 
over the integers and 

[ \ dt = TrJ2\f{h)f ■ 

Jm 

In mathematical terms, / is in the Paley-Wiener space PWi(R) consisting of all 
functions in L^(R) whose Fourier transform is supported in [—1,1]. More generally, 
in the theory of non-regular sampling a sequence A := {A}agA of points on the real 
line R is said to be sampling for PWi (R) if there exists a constant C such that the 
following sampling inequality holds 

2 [ \f{t)fdt<Y,\f{x)f<c[ \fit)\^dt 

^ Jr Jr 

for any any / £ PlUi(R), ensuring that the reconstruction of / is stable in the 
P^-sense. Corresponding results also hold in the higher dimensional setting where 
R is replaced with R" and the band [—1,1] with the unit-cube [—1,1]" (or more 
general any fixed convex body of volume one). By the seminal result of Landau 
[T7], a necessary condition for a set A to be sampling is that the corresponding 
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asymptotic density of points in R" (in the sense of Beurling) is at least equal to 
the Nyqvist rate I/tt", i.e. 


lim inf 
R —¥00 




> 



(uniformly over translations) for any smooth domain C K” assuming a uniform 
separation lower bound on the points in A. In one-dimension the reversed strict 
inequality is also a sufficient condition for sampling, but not in higher dimensions. 

By a rescaling, Landau’s density results may also be formulated in terms of the 
high frequency limit which appears when the frequency domain [—1,1]" is replaced 
with fc[—1,1]" for k large, i.e. PILi(R") is replaced with the corresponding Paley- 
Wiener space PVLfc(M”). In this context a sequence Ak := of sets of points 

on H." is said to be sampling for PII4(®.") if 


1 

C 



E 

A('=)eAfc 


/(AW) 


< c 


I 


\f{t)\^ dt 


for any / G PWk(RA) with the constant C independent of k. For the sake of 
simplicity if it is clear from the context we will omit the superindex k in 
and write simply A G A^. The corresponding necessary density condition on the 
sampling points may then be reformulated as 


lim inf 

R—^cx) 


#{Afc n n} 

fc" 


> 


-dt 


uniformly over translations for any domain C M” with |cin| = 0. (Landau’s 
setting corresponds to the case when Ak is of the form k~^A, but his arguments 
extend to this high-frequency setting). 

There is also a natural compact analogue of the Paley-Wiener setting on 
which is the one which is most relevant for the present paper, where R (or R") is 
replaced with the circle := R/27r. Then the role of PILfe(R) is played by the 
space Hk{S^) of all finite Fourier series on [0, 27r] with frequencies in [—fc, fc], i.e. the 
space of all trigonometric polynomials of degree at most k. A sampling sequence of 
finite sets of points Ak C in this setting is also called a Marcinkiewicz-Zygmund 
family [23) . In a similar setting to ours it has been studied in the sampling 
sequences for powers of a line bundle with positive curvature on compact complex 
manifolds. In this case there are precise estimates for the Bergman kernel so that 
Landau’s techniques carry through. 

From an abstract point of view the previous settings fit into a general Hilbert 
space framework where Hk{M) is a given sequence of Hilbert spaces of functions 
on a set M with reproducing kernels Kk{x,y). Then a sequence Ak of sets of 
points on M is said to be sampling for Hk{M) if the family of normalized functions 
ka := A"fc(-, A)/|| Arfe(-, A)|| for A G A^, form a frame in the Hilbert space Hk{M), in 
the sense of Duffin-Schaeffer [13], i.e.: 

^Il/f < E ia«A)P<C||/f , V/GHfe(M), 

Ag Afc 

which is equivalent to the sampling inequalities: 

( 1 . 1 ) h\ff< E '^f&Hk{M). 

^ tt. ^fc(A,A) 


where we will assume that C can be taken to be independent of k. 
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1.2. The present setting. The main aim of the present paper is to generalize 
the Landau type necessary density conditions for sampling on to a general 
setting where the Hilbert space Hk{M) consist of polynomials of degree at most 
k on an affine real algebraic variety M equipped with a weighted measure. We 
are not dealing with the very interesting problem of finding sufficient conditions 
for sampling multivariate polynomials. This and its numerical implementation is a 
very basic question in signal analysis, see for instance m and the references therein 
for the one-variable numerical sampling. 

Our setting is the following: by definition M is the variety cut out by a finite 
numbers of polynomials on K™ and Hk{M) is the space of polynomials of total 
degree at most k restricted to M and equipped with the norm 

Ibfc|li2(e-fc^ ) := / \pk\^ 

J M 

defined by a compactly supported measure fi on M and a continuous function (j) on 
M (referred to as the weight function). Following [4] we will refer to the pair (/i, (j)) 
as a ’’weighted measure”. In order that the latter norm be non-degenerate some 
regularity assumption has to be made on p. The affine case, i.e. when M = K™ is 
thus the classical setting for multivariate orthogonal polynomials. We will assume 
two regularity conditions: the Bernstein-Markov property and moderate growth, 
see Section lOI for the precise definitions. 

Our first main result in this general setting is 


Theorem 1. Let M be an affine real algebraic variety equipped with a non degen¬ 
erate measure p. and a weight function if. Assume that the pair (p, (p) satisfies 
the Bernstein-Markov property (EU and it is of moderate growth (12.21) . Then a 
necessary condition for a sequence of sets of points in M to be sampling for the 
space Hk{M) of polynomials of degree at most k, with respect to the weight kfi and 
measure p, is that 


( 1 . 2 ) 


hm mf — > peq 

fe->oo Nk ^ 

AGAfc 


in the weak topology on the measures on M, where peq denotes the normalized 
equilibrium measure of the weighted measure {p,4>) and Nk = dmi{Hk{M)). 


1.3. Sampling on compact real algebraic varieties equipped with a volume 
form. One disadvantage of the definition of the sampling inequalities (HU in this 
general setting is that it is of a rather abstract nature as it involves the reproducing 
kernel Kk{x,x) which in general is impossible to compute explicitly. On the other 
hand, only the asymptotic behaviour of Kk{x,x) as fc —>■ oo is needed and these 
asymptotics can often be estimated. Also, if peq is absolutely continuous with 
respect to Lebesgue measure than the condition (HU above may be written as 


(1.3) 


lim inf 

k—^oo 


#{Afc n 0} 


> 


peqiyi) 

peqiAI) 


for any smooth domain O.We will refer to the latter condition as the ’’pluripotential 
Nyqvist bound”. One particularly interesting case where Theorem H] can be made 
explicit is the following: 


Theorem 2. Let M be an n-dimensional affine real algebraic variety, which is 
non-singular and compact, let p be a volume form on M and let fi = 0- Then 
there exists a positive constant C such that the reproducing kernel for [Hk{M), p) 
satisfies 


(1.4) 


< Kk{x,x) < Ck'^ 

o 
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and thus {p, 4>) is non degenerate, it satisfies the Bernstein-Markov property and it 
is of moderate growth. Moreover, a necessary condition for a sequence Afc of sets 
of points on X to he sampling for Hk{M) is that the density of sampling points is 
at least equal to the density of the equilibrium measure fieq of M, as k ^ oo, i.e., 
the pluripotential Nyquist bound (IE3D holds. 

The definition of the equilibrium measure of M and more generally the equi¬ 
librium measure attached to a weighted measure will be recalled in Section 12.21 
As pointed out above this result thus generalizes the results in [53] concerning the 
case when M is the unit-circle. Moreover, the case when M is the unit-sphere 
corresponds to the case studied in [50] , where the signals in questions are spherical 
harmonics. While in all these special cases the equilibrium measure pLeq is explic¬ 
itly given by the Haar measure (since the corresponding Riemannian manifolds are 
homogeneous) the equilibrium measure of a general real affine algebraic variety 
appears to be of a highly non-explicit nature. Another generalization of the homo¬ 
geneous cases was considered in [53] , where the signals are “band-limited” sums of 
eigenfunctions of the Laplacian on a given compact Riemannian manifolds (M, g) 
and then the role of the equilibrium measure is played by the Riemannian volume 
form. 

It is not evident a priori that there are sampling sequences at all. This is assured 
with the following Bernstein type theorem; 

Theorem 3. Given a smooth compact real manifold M C R™ of dimension n, the 
following are equivalent: 

• M is algebraic 

• M satisfies a Bernstein inequality, i.e., for some q>l (or for all q> 1): 

l|Vtp||L,(M) < Cqdeg{p)\\p\\LpM)- 

• There is a uniformly separated Afc such that for some (all) q > 1 

f b|WM<4 E VpeiPfc(R™). 

Jm Jm 

This generalizes the main result of |9] where the case q = oo was considered. 

1.4. Sampling of multivariate real polynomials on convex domains. An¬ 
other instance where TheoremjTjcan be made more precise is the case were M = R”, 
/i is the Lebesgue measure restricted to a smooth bounded convex domain and 
(j) = 0. In this case the equilibrium measure is very well understood, see [T] and 
[Sj. It behaves roughly as dpeq — l/\/d(a;, dM)dV, (this will also follow from the 
asymptotics (11.511 below). 

Theorem 4. Let Ll be a smoothly bounded convex domain in R”. Then the repro¬ 
ducing kernel for {Hk(fil),dV) satisfies 

C A-" \ 

,I Vx G n. 

y/dix) J 

where d{x) denotes the distance of x G LI to the boundary of LI. Thus it satisfies 
the Bernstein-Markov property (EH) and it is of moderate growth (I23). Moreover, 
a necessary condition for the sequence Afc of sets of points on LI to be sampling 
for Hk{M) is that the density of sampling points is at least equal to the density of 
the equilibrium measure peq of LI, as k ^ oo, i.e. the pluripotential Nyquist bound 
(11.311 holds. 
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1.5. Interpolating sequences. A natural companion problem to that of sampling 
sequences are the interpolating sequences. In the same abstract point of view that 
we considered for sampling sequences we consider a sequence Hk{M) of Hilbert 
spaces of functions on a set M with reproducing kernels Kk{x,y) and instead of 
a frames we consider Riesz sequences of normalized reproducing kernels (see Sec¬ 
tion o for the precise definitions). 

Landau in m studied also these sequences in the Paley-Wiener space, and his 
observation was that locally if a sequence A is interpolating then its density should 
be smaller than the local density of the space. We can again use the ideas inspired 
in [22] to deal with the case of polynomials in real algebraic varieties. 

Our main result to this problem is 


Theorem 5. Let M be an affine real algebraic variety equipped with the Lebesgue 
measure. Then a necessary condition for a sequence Aj, of points on M to be 
interpolating for the space Hk{M) of polynomials of degree at most k, is that 


( 1 . 6 ) 


lim sup —— 

fc^OO Nk 


^ yeq 


in the weak topology on M, where y,eq denotes the normalized equilibrium measure 
of M and Nk = dim{Hk{M)), i.e. the following reversed pluripotential Nyqvist 
bound holds: 


lim sup 

k—¥C!0 


#{Afc n 0} 




< 


for any given smooth domain Q in M. 


yeq (O) 

k-eq{M) 


1.6. Discussion of the proof of Theorem|2j Let us make some brief comments 
on the circle of ideas involved in the proof of the previous theorems. First of all, since 
the sampling points uniquely determine a polynomial pk on M the total number 
#Afc of sampling points at level k is of course at least equal to the dimension Nk of 
Hk{M). As a well-known guiding principle the necessary conditions for sampling 
should come from a localized version of this argument saying the asymptotic lower 
density of sampling points should at least be given by the “local dimension” of 
the Hilbert space were p-k '■= e~^^dy, i.e. by the leading 

asymptotics of times the function 


Bk{x] 


Nk 

-E 


(k) f \ 

\Pi (x) 


i=l 


^-k4>{x) 


where {pPE)} is any orthonormal base in the Hilbert (note 

that integrating Bk{x) with respect to dV indeed gives the dimension Nk of Hk{M)). 
The independence of the choice of base follows from the following extremal repre¬ 
sentation of Bk{x): 


(1.7) 


Bk{x) 


sup 

Pk&Hk{M) 


|pfc(3:)|^e-^'^W 

SM\Pk? 


From the point of view of general Hilbert space theory Bk{x) may be written as 
Bk{x) = Kk{x,x) where Kk{x,y) is the reproducing kernel of the Hilbert space 
{Hk{M), |Ml 2 , 2 (d/ife))’ i-®- i-ii® kernel of the orthogonal projection from L^{dy,k) to 
Hk{M). By the general results in |3] 


( 1 . 8 ) 


N^. ^Bk{x)dy,k yieq/PeqiM) 


weakly on M as fe —>■ oo, which in the view of the guiding principle above thus gives 
a strong motivation for the previous theorems. However, this guiding principle does 
not seem to hold in all generality and it has to be complemented with some further 
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asymptotic information of the full reproducing kernel Kk{x,y). This is already 
clear from Landau’s classical proof in the Paley-Wiener setting on K- HZ], where 
the decay asymptotics of Kk(x^y), away from the diagonal are needed in order to 
construct functions fk which are well-localized on a given domain 17 (using suitable 
Toeplitz operators). Moreover, Landau’s approach also relies on certain submean 
inequalities for fk which pose difficulties in our general setting. Instead we use a 
new approach to proving necessary conditions for sampling, which is inspired by 
and m, where we reduce the problem to establishing two asymptotic properties 
of Kk{x,y) (given the convergence of the Bergman function to the equilibrium 
measure): 

• A growth property of Bk 

• A weak decay property of \Kk{x,y)\ away from the diagonal 

The interpolation theorem [S] is proved in a similar way, but by replacing the 
moderate growth property with a Bernstein type inequality, see Theorem [HI 

One interesting feature is that although the statement of the problems studied are 
purely real, all the proofs rely on the process of complexification with one important 
exception: the off-diagonal estimate on the Bergman kernel (Theorem 1171) which 
exploits, in an essential way, the real structure. 

1.7. Further relations to previous results. As explained above one important 
ingredient in the proof of Theorem [T] is the asymptotics for Bk{x) in formula (HU) 
established in [1], which holds generally under the Bernstein-Markov assumption 
on (/r, (jj). In turn, the latter result can be seen as a consequence of a very general 
result in H] giving the convergence towards the equilibrium measure of (^, </>) for the 
normalized Dirac measure associated to a sequence of Nk points under the condition 
that the points are asymptotic Fekete points for the weighted set (Supp(/i),(()). It 
may thus be tempting to try to deduce Theorem [T] in the present paper directly 
from the general convergence results in [1] by the following tentative procedure: 
one removes points from a sampling sequence Afc until one arrives at Nk points, 
while keeping the sampling property. But as shown by a counter example in ESI 
Example 2] such a procedure is doomed to fail, already in the homogenouous case 
of the circle. 

The point-wise asymptotics for Bk{x) in TheoremE|(formula (11.41) 1 can be seen as 
an improvement - in the special case of a real algebraic manifold - of a classical result 
for regular compact subsets complex space going back to Siciak and Zaharyuta 
giving that i3fe(a;)^/^ —>■ 1 point-wise on M (this latter classical result was given a 
i9—proof by Demailly, m which can seen as a precursor to our proof). In view of 
the weak asymptotics (11.81) and the bounds in formula (12.31) below, it seems natural 
to conjecture that Bk{x)/k^ in fact converges pointwise (in the almost everywhere 
sense) to the density of the equilibrium measure of M (and similarly in the 
setting of a convex domain; as in the one dimensional setting, see EHI). This would 
be a real analog of the point-wise asymptotics for Bk {x) in E| i where the role of M is 
played by a complex projective variety endowed with a hermitian holomorphic line 
bundle (the case of positive curvature is a fundamental result in complex geometry, 
due to Bouche [7] and Tian EH). 

On the other hand, in the line bundle setting, the asymptotics for \Kk{x,y)\^ in 
Theorem HZI are only known in the case of a line bundle with positive curvature, 
but it seems natural to expect that they hold in general (see El for some results 
in this direction in connection to the study of fluctuations of linear statistics of 
determinantal point processes). 

Finally, we recall that in the one-dimensional case there is a vast litterature on 
various asymptotic results for orthogonal polynomials, in particular in connection 
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to random matrix theory. For example, the asymptotics of the scaled reproducing 
kernel k~^Kk{x + ^,x + ^) have been established, in connection to the question 
of universality, under very general condition on a given measure /r on the real line, 
when X is a fixed point in the “bulk” of the support of /r; see for example [19] (and 
similar scaling result holds at the “edge” of the support). However, there seem to be 
very few results in the higher dimensional setting (but see [29] and [16] , where the 
case of the ball and the simplex in is settled). It would be interesting to extend the 
asymptotics in the present paper to study similar universality questions in higher 
dimensions and we leave this as a challenging problem for the future. 

Acknowldegements: The authors are grateful to Ahmed Zeriahi for very useful 
discussions regarding the Bernstein inequality. 

2. PlURIPOTENTIAL theory and ASYMPTOTICS OF REAL ORTHOGONAL 

POLYNOMIALS 

2.1. Setup. Let M be an n-dimensional affine real algebraic variety, which is non¬ 
singular and compact. In particular, M is the common zero-locus of a collection 
of real polynomials pi,...,pr in R™. We denote by Hk{M) the real vector space 
consisting of the functions on M which are restrictions of real polynomials in R™ 
of total degree at most k. We will also consider the “complexifications” X and 
Hk{X) of the real variety M and the real vector space Hk{M), respectively. More 
precisely, X is the complex algebraic variety in C™ defined by the common complex 
zeros of the ideal defining M and Hk{X) is the complex vector space consisting of 
restrictions to X of polynomials in C™ of total degree at most k. Then M is 
indeed the real part of X in the sense that it consists of all points in z in A such 
that z = z and real vector space Hk{M) is the the real part of Hk{X) in the 
sense that it consists of all pk in Hk{X) such that Pk = Pk (restricted to M). 
Denoting by X the closure of X in P™, which defines a compact (possibly singular) 
complex projective algebraic variety we will, in the usual way, identify Hk{X) with 
the space H°{X,Ox{^)^^) of all global holomorphic sections of the line bundle 
Ojf(I)®^ ^ A. 

We will denote by Kk{x, y) the Bergman reproducing kernel of Hk equipped with 
the L^-norm induced by a given weighed measure (/i, (j)) and moreover we will use 
the notation Bk{x) = Kk{x,x)e~^^''^\ 

Example 6. Let M be the unit-circle realized as the zero-set in R^ of p{x,y) = 
— I. Setting x = cos0 for 0 G [0, 27r] we may identify Hk{M) with the space 
iLfc([0,27r]) of all Fourier series on [0, 27r] “band-limited” to [—k,k] i.e. spanned 
by 1 and cosmO and sinm0 for m G [l,fc] fl Z. More precisely, iLfc([0, 27r]) is the 
pull-back of Hk{M) under the corresponding map from [0, 27r] to M. To see the 
relation to the more standard setting where M corresponds to the unit-circle in 
C with complex coordinate r (equal to e*® on 5^) we note that the embedding F 
of C* in given by 

z := {t + t~^)/2 w = {t — T~^)l2i 

maps C* to the complex quadric A cut out by p(z, w) := z^ + w'^ — 1 and the 
unit-circle in C* is mapped to the real part M of A. Indeed, r G 5'^ if and only 
if r = T~^ iff z(t) = 5Jt(= cosO) and w{t) = 3t(= sin0). The pull-back of Hk{X) 
under F is the space of Laurent polynomials on C spanned by the monomials r™ 
for m G [—fc, k\ fl Z. Hence, the real part of F*Hk{M) is indeed spanned by I and 
cosm9{= 5ftr™) and sinm0(= 3 t™) for m G [I, k] fl Z. 

Definition 1 (Bernstein-Markov). The standard assumption on the pair (/r, (p) is 
that it satisfies the Bernstein-Markov property (with respect to the support of p ), 
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which may be formulated as the property that the reproducing kernel Kk{x, y) have 
sub-exponential growth on the diagonal, i.e. for any e > 0 there exists a positive 
constant Ce such that 


Bk{x) < 


( 2 . 1 ) 


uniformly on the support of y,. 

For our general results to hold we will also need the following technical regularity 
of the growth assumption on the reproducing kernel along the diagonal: 

Definition 2 (Moderate growth). We say that Hj. has a reproducing kernel with 
moderate growth if 


Kk+i{x,x) < CKk{x,x) 


( 2 . 2 ) 


on the support of /i. More precisely, for our purposes the constant C may be 
replaced by any sequence with growth of the order o(fc). 

Anyway, in all examples that we are aware of the constant C in (12.2|) may actually 
be replaced by a sequence tending to one as fc —>■ oo (which, by iteration, actually 
implies the Bernstein-Markov property (EU). All the measures y supported in M 
that we will consider are non-degenerate, in the sense that ||/|| ^^OifO^/G 
Hk{M). Otherwise if the support of y is contained in the zeros of a non vanishing 
polynomial in Hk{M) one may replace M with a subvariety. 

While Definition [1] is standard (see [ 3 ] and references therein), Definition [ 5 ] ap¬ 
pears to be new. We expect it to hold in great generality and we will establish it 
in the situations relevant to the present paper. 

2.2. The extremal function attached to a real affine variety. Recall that 
the Lelong class £(C"^) is the convex space of all plurisubharmonic (psh, for short) 
functions (/> on C'” with logarithmic growth, in the sense that (/> < log(l -I- |zp) -hC. 
The restriction of this space to X will be denoted by £(X) and it may be identified 
with the space of all (singular) metrics on the line bundle 0(l)x —t X with positive 
curvature current, see [ 3 ] and references therein for further background. The Siciak 
extremal function (sometimes called the equilibrium potential) of a compact and 
non-pluripolar subset K oi X and a weight if G C(IC) is the function in £(X) 
defined as the upper semi-continuous regularization of the envelope 


'VK,i>{x) = sup {4){x) : (j) < Ip on K} 

4>G/:{x) 


and the weighted set {K,ip) is called regular if = VK,i)- The Monge-Ampere 
measure 


yK,i) '■= MA{vk,iIi) 


is called the (pluripotential) equilibrium measure of {K,ip) and it is supported on 
K (when there is no risk of confusion we will write ■= fJ-eq) 

In the following we will take K := M as above, which is thus embedded in the 
complex affine variety X in such a way that M = A fl {j/ = 0}, where y denotes 
the imaginary part of z G C". We will also take ip = 0. 

Proposition 7. Let M be an n-dimensional affine real algebraic variety, which is 
non-singular and compact and denote by vm its extremal function, defined on the 
complexification X of M. Then M is non-pluripolar and regular. Moreover, there 
exists a constant C such that 


(2.3) 
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in a neighborhood of M in X. In particular, the equilibrium measure piM is abso¬ 
lutely continuous with respect to the Lebesgue measure dVu on M and its density 
is bounded from above and below by positive constants: 

(2.4) j^dVM < g-M < DdVM 


on M. 


Proof. The lower bound on vk follows from a simple max construction. Indeed, we 
may after a scaling assume, as before, that |z| < 1 on M and then set 


^ := 


max{|?/|/C', log|zp} 
\og\z\'^ 


if 1^1 < 2 
if l^l > 2. 


with C sufficiently large to ensure that "0 is continuous along \z\ = 2. Since, 
'll: = \y\/C close to M the function i/i is a contender for the sup defining vm and 
hence if < vm, which proves the lower bound in (12.311 . The proof of the upper 
bound is more involved: We know that vm G L^^{X) because M is algebraic, see 
|26) . Moreover there is a distance d such that g{r) = sup^^x:d{z,M)=r'^M(z:) is 
convex in r in [0, ^], see the proof of Theorem [TU] for details. Therefore g{r) < Cr 
and the upper bound in (12.31) follows. 

Anyway, for the proof of the lower bound in (Ha we will only need the lower 
bound in (12.31) . Note also that combining (d and the asymptotics dUD imme¬ 
diately gives the inequalities (12.41) (which also follow from the inequalities (I2.3|l by 
the comparison principle for the Monge-Ampere measure, see [TJ Lemma 2.1]. □ 


2.3. The proof of the lower bound on Bk in (11.41) . Denote by BkvM the 
Bergman function on X defined by the L^-norm on Hk{X) induced by the measure 
The idea of the proof is to first show that 

(i) BkvM > on M, 

and then that 


(ii) 


Bk > Ck-^Bky^ on M. 


This would clearly imply the result in question. However, for technical reasons we 
will only show a slightly weaker version of these inequalities (needed for Q) where 
Vm is replaced by 


Vm ■= wm(1 - e)-\-eil) 

where if is a, continuous function in C{X) such that if = v^^jC in a neighborhood 
of M. Here e is a sufficiently small positive number which is fixed once and for all. 
To see that such a function if exists we may after scaling assume that jzj < 1 on M 
and then simply set if := maxju^/C, log \z\'^} when jzj < 2 and (f := log \z\^ when 
l^l > 2. The constant C is taken sufficiently large to ensure that if is continuous 
at | 2 ;| = 2. 

Let us start with the proof of ((il). To simplify the notation we will assume that 
n = 1 (but the general proof is essentially the same). To this end fix a point in 
M and introduce local holomorphic coordinates z on U in A, centered a the fixed 
point, such that M = {y = 0} locally, i.e. on U (not to be confused with the global 
coordinates on C™ and M", respectively). The idea is to first construct a local 
function fk, holomorphic on U such that 


(2.5) 


IA-(0)|^ 


lulfkr^^dV 


> k^/C 


and then perturb fk slightly to become a polynomial pk by solving a global d- 
equation on X with an L^-estimate. 
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There is no loss of generality assuming that /fc(0) = 1. Working in a local 
coordinates and reescaling (1^ it is enough to prove that there is a function / G 
'H(C) such that /(O) = 1 and 

Jc 

then fk{z) = f[kz) will satisfy (12.51) . The function f{z) = smc'^(Cz/2) has the 
desired properties. 


2.3.1. Modification and globalization. Let now y be a smooth cut-off function sup¬ 
ported on U (say equal to one on C//2). In view of standard globalization arguments 
the problem with the present setting is that d of the global function xfk on X does 
not have a small weighted L^- norm (compared to the weighted norm of fk). The 
reason is that is only well localized in the y-direction. To bypass this dif¬ 

ficulty we will instead replace vm with v%i and modify fk accordingly as follows. 
First observe that, by definition, 




|2 -kvM{i-e) p-kei\y\^ 


for any gk- We next observe that 
gives 


e = e and hence setting gk := 


and yfc(O) = /fc(0). In particular, 


2 -kie\v\ _\f\2 -2k<i\z\ 


= IM <\fk\^ 


|gfc(0)| 


- '"n ' '• - > - 


i/fe(o)r 


I /■ |2 


> fcVa 


Here the optimal constant is slightly smaller than the previous optimal C, but 
on the other hand we have gained a Gaussian factor that we will next exploit. The 
point is that dfxgk) = dx9k is supported where j^l > 1/4 and hence 


\d{x9k) 


12 ^ — kv 


^dV < C 


I 


l/4<p|<2 




C 


I 


l/4<p|<2 




Estimating the Gaussian factor e -^ith its sup, i.e with e ^2e/4 gives 

the bound 


[ \dix9k)\ 

Ju 




0(6-^'=) f \fk\ 

Ju 


2 g-fewMll-c) 


Here and henceforth 0{e~^^) denotes a term which is exponentially small in k 
(recall that that e is a small number which is fixed once and for all). 

With this local estimate in place we can now apply a standard globalization 
argument: using L^-estimates for d on the line bundle 0(l)x over X, or more 
precisely (if the latter variety is singular) on its pull-back to a smooth resolution of 
X) there exists a smooth function Uk on such that pk ■= 9k — Uk is in Hk{X) and 


duk = d{xgk), [ \ukf e [ \d{x9k)\^ e 

J X Ju 

(strictly speaking to apply L^-estimates we have to slightly modify the weight v\j 
with a fc-independent term to ensure that the corresponding metric on the line 
bundle kO{\)x has a sufficiently large uniform lower bound on its curvature form, 
but this only changes the L^-estimates with an overall multiplicative constant. 
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which is harmless). This is a standard procedure; for a precise statement which 
also applies in the singular setting see, for example [31 Section 2]. 

By the previous estimate this means that 

Jx Ju 

Moreover, applying the mean value property for holomorphic functions on a small 
coordinate ball then gives 


Uk{0)<Ck^ [ \ukfe-'^^‘^dV <0{e-^'^) [ |/fc|" 

Jx Ju 

Hence, 

bfc(o)l^ ^ |gfc(o) - ^ 

Jx\Pk\'^e-'=^MdV Jx\xgk-Ukf e-'^'^MdV ~ 

lu Ixflfcl^ e~^'"MdV + 0(e-'5'=) |/fc|^ 


But |gfe(0)|^ = |/fc(0)|^ and \gk\'^ e < \fk\^ e Moreover, as explained 

above \ fk\^ _ 0{k^) and hence we get just as above 


\Pk{0)? 




which concludes the proof of the bound ([il) 


Bkvi, > Ck^^ on M. 


2.3.2. The inequality between Bk and BkvM ■ First observe that it is enough to prove 
the following lemma where now y denotes the imaginary part of z £ C" (so that 
Xn{y = 0} = M}) : 

Lemma 8. Let Uk be the set of all points in X such that |?/| < 1/fc (which defines 
a neighborhood of M in X). Then there exists a constant C such that 

for any polynomial of total degree at most k. 

Indeed, since the function vm on X is comparable to |y| close to M (by Theo¬ 
rem El) and in particular kv\j is uniformly bounded on Uk-, we then get that 

[ \Pk?dV<C'^[ \pk?e-’^^MdV 
J M ^ Jx 

It follows immediately that 

Bk > k^lC'Bk.^^ 

on M, which combined with the inequality ([il) thus concludes the proof of the lower 
bound in (1121), given Lemma [51 to whose proof we next turn. 

For any x £ M there are constants C and rp such that for any r < tq, 

|/(x)|^<4r/ \f{v)\''dVx{y). 

^ J XnB(x,r) 
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for any / holomorphic in X. In particular if we integrate over x € X a polynomial 
of degree k taking r = 1/k we get 


f \Pk?dVM<Ck'^'^ f \pk{y)?VM{B{y^/k)f^M)dVx{y) 

M JUk 

< Ck^ [ \pk\^dVx < [ \Pk\^dVx. □ 

Juk Vol{Uk) Juk 


2.4. The L‘^ Bernstein inequality. Let M be a smooth compact algebraic variety 
in R"* of dimension n. 

Given a polynomial p € VkiW^) and x G M we denote by Vtp{x) the tangential 
gradient of p along the manifold M. The following Bernstein type inequality holds: 


Theorem 9. Let q G [l,oo], then there is a constant Cq such that 

l|Vtp||L,(M) < Cqdeg{p)\\p\\Li(M)- 

The case q = oo was proved in [Q. We prove now the case <7 = 1 and the others 
follow by interpolation. 

Let X be a complexification of M, i.e, an algebraic variety in C'" such that 
M = X (1 R'". We denote hy Ur C X the neighborhood of M defined as Ur = {x G 
X : d{x, M) < r}. By the Cauchy inequalities we have that for any x G M and any 

feuix): 

\^tf{x)\ < f \f{y)\dVxiy). 

X JB{x,r) 

and integrating over M we have 

[ \XJ{x)\dVM < [ [ \fiy)\dVxiy)dVMix) < ^ / \fiy)\dVx{y). 

Jm X J M J B(x,r) X Jjj^ 

Therefore Theorem [5] follows from the following result: 


Theorem 10. There is C > 0 such that for all polynomials pk of degree k, the 
following inequality holds: 



\Pk\dVx < Ck-'^ [ \pk\dVM. 
Jm 


Proof. In order to estimate the integral over Ui/k we will integrate along surfaces 
surrounding M. These surfaces will be level sets of plurisubharmonic functions 
with Monge-Ampere 0. In this setting there is a generalization of Hadamard three 
circles theorem due to Demailly that will be used, see [T2]. In order to use this 
technique we need that the psh-function that defines the level sets is smooth out of 
M, and its square must be smooth. We can use the function provided by Guillemin 
and Stenzel in |15j in their study of Grauert tubular neighbourhoods around real 
analytic manifolds. We present the setting: 

Take if a plurisubharmonic function in a neighborhood U of M in Ai defined as 

if{z) = d{z,M), 


where the distance d is given by a metric provided in a Grauert tubular neighbor¬ 
hood U as in m\- The function satisfies {dd'^tjj)^ = 0 in U \ M, is a real 
analytic Kahler potential in U and {dd'^{tjj‘^))^ is comparable to the volume form in 
AT in a neighborhood of M. 

We use the same notation as in m- Consider the pseudospheres S'(r) = {z G 
U; if{z) = r, r > 0 and the positive measures pr supported on Sr that are defined 
as 

f h{dd‘^f,Y-^ Adfif 

Js{r) 
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for any h € C{U). When r = 0, then we define 

/2o(/i) := f h{d(fij))‘^. 

J M 

We have that ^rih) is continuous for r > 0 with —?► fio{h) as r —>• 0+ see [H 

Theorem 3.2]. 

Moreover {dd^ip)'^ which is supported on M is comparable to the volume form 
in M. This is so because at any point z € M we can write local holomorphic 
coordinates such that M corresponds to z G C" : = 0. In this coordinates '4>{z) 

is comparable to jSzj and since is the Lebesgue measure on M" then by 

the comparison principle for the Monge-Ampere measure, see [TJ Lemma 2.1], the 
measure /io is locally comparable to the volume form and M being compact it is 
globally comparable. 

Take the psh function V = loglpfe], then [T^ Corollary 6.6(a)] says that the 
function 

u(r) = logr > 0, u(0) = log p,o(e^) 

is convex and increasing in r. We fix i? > 0 such that Sr belongs to the tubular 
neighborhood U. The convexity of u implies that for any r > 0 

(2.6) u{r) < + u{R)^. 

We have that 

u{0) = log \pk\dp,o. 

Jm 

We are going fix R such and estimate u{R). Since pk is a polynomial of degree k 
we have that by the Bernstein-Walsh estimate 

suplpfe] < supsup Ip], 

Sr Sr M 

where (j^M is the Siciak extremal function defined as 

= max|o,sup{^j^^^loglp(z)l : pG ^^(0™),deg(p) > 0,suplpl < 1}| 

It is a well-known theorem of Sadullaev, see m that if X is algebraic then 
4>k G for any non-pluripolar compact set K relative to X. Certainly M is 

non-pluripolar relative to X since it is totally real. Therefore: 

suplpfc] < C'^sup 1 p1 
Sr m 

Moreover, we need the following Bernstein-Markov type inequality: 

(2.7) sup\p\<Clj f \p\. 

M Jm 

This is easier than the standard Bernstein-Markov property since we are not re¬ 
quiring that Cm is close to 1. In our case (EZD is a special case of [51 Theorem 4.1]. 
Finally supg^ Ipfe] < 1 p 1, and we have that u{R) < (C -I- log ||p||Li(M))fc- 

Therefore if r = 1/A: and using the convexity (12.6|) we deduce that 

uiyl/k) < u(0) -I- C, 

Since u{r) is increasing we have that for any r <\/k 


/ \Pk\dpLr <C 


IPfeM^o- 




14 


ROBERT J. BERMAN AND JOAQUIM ORTEGA-CERDA 


But the measures fir desintegrate the form ^ ^d'ip Ad‘^'tlj, see [T2J Propo¬ 

sition 3.9], and we have that 


rl/k p p 

/ \pk\dfjir=l A di/’A 

Jq J S{r) Jipcil/k 


Moreover 

i/;(2)’^-^(dd=V')””^ Adi) A d^ii = (dd=(i/;2))" 

But "0^ is a real analytic Kahler potential in X, see [15]. Thus {dd‘^{ip'^))'^ is 
equivalent to the original volume form Vx in a neighborhood of M in X: 


nl/fc 


' '0<1//C 


\Pk\dVx ^ [ r'^ ^ f \pk\dfir 

Jo Js{r) 


Remark. With the same proof, for any 1 < g < oo, 

[ bfclWjf <fc-” / bfejWM. 
ic/i/fc Jm 


It is also true that 

sup jpfej < sup jpfcj. 

Ui/k M 

The proof is the same, but instead of m Corollary 6.6(a)] one uses that 

u(r) = suploglpfe], 

Sr 


is a convex function of r, see [TU Corollary 6.6(b)]. 


□ 


Remark. For any x £ XI we consider a ball Bx{x, 1/k) in the complexified manifold 
X. By the submean value property we have that 

\Pk{x)\'^ < [ bfc(y)P dVx{y) < /" bfcP dl4f 

JBx{x,l/k) Jm 

Therefore Kk{x,x) < fc" and we have proved the upper inequality in (11.41) . We 
include the argument for completeness but this upper bound is well known and it 
follows from the arguments in m- 

It is also possible to prove a converse result to Theorem [HI 


Theorem 11. Let M be a smooth compact submanifold in K™. If there is a con¬ 
stant C > 0 such that for some g G [1, oo], 

l|VtpllL9(M) < C'deg(p)l]pl]i,(M), 
for all polynomials p G VfMX'), then M is algebraic. 

Proof. We will need a definition 

Definition 3. A sequence of finite sets {A^} is an e-net if it is uniformly separated 
and 1 < Z^AGAfe XB(\,e/k){x) < C, for all x € M and k>0 

By an application of the Vitali covering lemma it is possible to construct e-nets 
for arbitrarily small e where the constant C = Cm depends on the dimension of M 
but not on e. 

Given an e-net A = A(e) we denote by Ik = #Afe. We may define: Tk ■ Vk{M) —>• 
as 

Tk{p){A) = PBuXe/k) '■= T pdVM VA G Afe. 

JBM{\,e/k) 
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We will prove now that if £ is small enough then 


( 2 . 8 ) 


f i»i' < ^ E iA(p)(A)r 

JM « 


AGAfc 


If this is the case, then Tk is one to one and dim(7^fe(M)) < Ik — where n = 
dim(M) which is much smaller than fc"*. Thus M is algebraic. 

Let us prove (12.811 . 

[ < -^ X! +X! / \pix) -pBM{\,e/k)\‘‘dVM- 

Jm ^ ^^A jBM(X,e/k) 


AeAfc 

By the Poincare inequality: 


[ \pix)-PBM{Ke/k)\'^dVM f IVtpl'^dVM- 

J BM(\,e/k) J BM(^,e/k) 


' Bm 

By Theorem [5] 


f M 


IPkl^^dVM < T E \Tkip)iXW+e [ \pk\‘^dVM. 

* aga. Jm 


and if e is small enough then (12.811 follows. 


□ 


2.5. Applications. We will use now the Bernstein inequality to get some more 
information on sampling and interpolation sequences of finite sets. 

Definition 4. A sequence of measures {pk}k is said to be a uniformly sequence of 
Carleson measures if there is a C > 0 such that 

(2.9) [ \pk\‘^ dfik < C\\pk\\‘^ , 

Jm 

In M we consider the balls defined by any Riemannian metric. 

Proposition 12. A sequence {p,k}k is a uniformly sequence of Carleson measures 
if and only if there is a C > 0 such that 

Pk{B{x,l/k)) < Cjk^ for allxGM,kG N. 

Proof. The necessity follows from testing (1^ against normalized reproducing ker¬ 
nels. For any x € M, let Kx,k{y) = Kk{x, y)/y/Kk{x, x). Then it is clear that for 
all y G M: 

\Kx,k{y)\ < \{Kk{x,-),Kk{y,-))\/\/Kk{x,x) < ^/Kk(y,y) ~ 

and \Kx,k{x)\ = \jKk(x,x) ~ On the other hand since M is algebraic we have 
the classical Bernstein inequality, see [5]: 

sup ||VtA£ai,fc|| < A: sup^ . 

M M 

This means that there is a (5 > 0 such that 

Wx,k{y)\ > K,x,k{x) ~ Vy G Bm{x, S/k). 

Therefore if we test (1^ with Kx,k we get that y.k{BM{x,S/k)) < I/fc". 

On the other direction for any x G M we consider a ball Bx{x, r) the ball in the 
complexified manifold X. By the submean property of the holomorphic functions 
/ G 'H(A), we have that for r < ro, \f{x)\ < fBj,(x,r) \fiy)\dVx[y)- Thus 

\Pk{x)\'^<i \pk{y)\'^dVx{y)^k^^ [ 


\Pk{y)? dVx{y). 
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Finally, thanks to Theorem [TUI 

[ \pk{x)\'^ d^j.k{x) < [ /" \Pk{y)\'^ dVx{y),dpk{x) < 

Jx JxeX JBx{x,l/k) 

k^^\Pk{y)\^MBx{y, i/k) n x)dVx{y) < 

k'"\Pk{y)\'^dVx{y) < [ \pk{x)\'^ dVM{x). 
Jx 


I 

L 


yeX,d{y,M)<l/k 

f 

yeX,d{y,XI)<l/k 


□ 


An immediate corollary is the description of the sequences that satisfy the left 
hand side inequality of the sampling sequences, that is a Plancherel-Polya type 
inequality. 

We say that a sequence of finite sets is uniformly separated if and only if 
there is an e > 0 such that d(A, A') > e/k for all A ^ A', A, A' G A^. 


Corollary 13 (Plancherel-Polya type inequality). The sequence of finite sets Ak 
is a finite union of uniformly separated sequences if and only if there is a constant 
C > 0 such that 


1 

k'^ 


E IpWI' ^ C' 



Vp G Vk- 


Proof. Apply Proposition [T^] to the measures Pk = -^ SaeA^ This implies that 
the Plancherel-Polya type inequality holds if and only if ff{Ak fl B{x, 1/fc)} < C 
uniformly in x and k. That is, the sequence A^ is a finite union of uniformly 
separated sequences of sets. □ 


Once we have a Bernstein type inequality the following Proposition is standard, 
see EH Proposition 5, p. 47] and it allows to reduce our considerations to uniformly 
separated sequences. 

Proposition 14. If Ak is a sampling sequence then there is a uniformly separated 
sequence of subsets A], C A^ such that A], is still a sampling sequence. 

It is also completely standard that 

Proposition 15. If Ak is an interpolating sequence then is is uniformly separated. 

Theorem [TOj can also be used to provide a sufficient condition that assures the 
existence of sampling sequences. More precisely, we say that the sequence Ak is an 
e-net if it is uniformly separated and 1 < X^AeAfc XB{x,e/k){x) < Cm, for all x G M 
and k > 0, the constant Cm depends on M but not on e. By an application of the 
Vitali covering lemma it is possible to construct e-nets for arbitrarily small e. 

Proposition 16. There is an eg such that any sequence Ak that is an e-net with 
e < eg is a sampling sequence. 

Proof. Take one e-net Ak = Afc(e). Then 

[ \Pk\dVM< E f \Pk\dVM< 

JM AgAfc dBM{X,e/k) 

< E bfe(^)ll-®M(A,e/fc)| + E [ \Pkix) - PkW\dVM < 

AeAfc AGAfc BM{\e/k) 

c-n _ p 

^ E k;i\P'^W\+ E ^\ytPk{<:x)\\BM{X,e/k)\, 

A^Afc A^Afc 
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where Ca £ BMiX,e/k) is such that |VtPfe(CA)| = swp^^BM{\,e/k) \^tPkix)\. By the 
Cauchy inequality, if we take a ball Bx{X, 1/fc) in the complexification X of M, we 
have 


|Vip(CA)| < fc / \pk\dVx. 

JBx{X,l/k) 


Therefore, 


[ \Pk\dVM < £ E + E [ \pk\ dvx. 

JBx{X.l/k) 

Since is an er-net there are at most Ce"" points of Ak in any given ball Bx (x, 1/fe) 
of center x € Ui/k- Thus, 

[ \Pk\dVM<^ E bfe(A)|+efc" f \pk\dVx. 

We use now Theorem [TO] to control the right hand side integral. If we take e small 
enough we can absorb the integral in the left hand side and we get 

[ \Pk\dVM < ^ \pk{X)\ 

* aga. 

The L°° version: supj\,f \pk \ ^ sup^^ \pk \ follows immediately by the Bernstein 
inequality proved in [^ , if £ is small enough. By interpolation we get that for any 
q G [l,oo) 


A^Afc 


The reverse inequality 


^ E ^ [ IPkl'^dVM, 

\ /- A . M 


Ae Afc 


follows from Corollary 1131 since is uniformly separated. 

We can finish now the proof of Theorem |2| 


□ 


Proof. We have already proved that the algebracity of M is equivalent to the Bern¬ 
stein inequality, this is Theorem [HI and [TTJ Moreover Proposition [TO| proves that 
compact algebraic manifolds have uniformly separated sampling sequences. So we 
only need to check that if there are such sequences then M is algebraic. This is 
proved in a similar way to Theorem |TT| We denote by Ik = #Afc as before. Define: 
Rk : Vk{M) —>• as 

Rk{p){X) = p{X) VAeAfc. 

Clearly, since we have the sampling property, Rk is one-to-one. Therefore dim(Pfc(Af)) 
Ik- Moreover since Afc is uniformly separated, then Ik < fc". This implies that M 
is algebraic. □ 

2.6. A general ofT-diagonal estimate on the reproducing kernel. 

Theorem 17. Let M be an n-dimensional affine real algebraic variety (possibly 
singular), pk a sequence of non-degenerate finite measures on M with support con¬ 
tained in a compact of M and denote by Kk{x,y) the reproducing kernel for the 
space Hk{M), viewed as a subspace of L^{M,pk)- Then there exists a positive 
constant C such that 

f I 2 2 

/ dfik{x)<S>dfikiy)\x-y\ <C/k 

JMxM ^ 
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Remark. Observe that if we pick /i^ = e the theorem covers the weighted 
setting as-well. 

Proof. Given a bounded function f on M we denote by Ty be the Toeplitz operator 
on Hk{M) n L'^{M,fik) with symbol /, i.e. Tf := 11^ o /• where 11^ denotes the 
orthogonal projection from LF‘{M,^k) to Hk{M), i.e. Tf is the Hermitian operator 
on Hk{M) determined by 

i'^fPkiPk)~ 

for any pk G Hk{M). The following is essentially a well-known formula 

TrTf -TrTf 2 = ^ f \Kkix,y)\'^ dpk{x)'S) dpk{y) (fix) - f{y)f 

^ Jmxm 

We provide nevertheless a proof for convenience of the reader: 

Claim. Let H be a reproducing kernel Hilbert space with kernel K, then for any 
bounded symbol f we have 

II \f{^) - f{y)\^\K{x,y)\^ = Tr(2T|/|2 - Tf o Tf - Tf o Tf). 

Proof K{x,y) = T,ufnix)My) 

Tf{g){x)= ! K{x,y)f{y)g(y). 

We compute the traces of T|j |2 and of Tf o Tf. 

Tr(r|/|2) = ^(/„,r|^|2(/„)) = 

n 

f / K{x,y)\f\^{y)fn{y) = [[\K{x,y)\^\f\^{y). 

^ J X J y JJ 

Thus 

Tr(%| 2 )= II \K{x,y)\^\f{x)\^ = J \Kix,y)W{y)\^. 

Now 

Tr{Tf oTf) = Y^ f fn{x) f K{x,y)f{y)Tf{fn){y) = 

n •'V 

= [ K{y,w)f{w)fn{w) = 

= 111 K{x,w)K{y,w)K{y,w)f{y)f{w) = II \K{y,w)\'^f{y)f{w). 

Similarly 

TriTf-oTf) = II \K{y,w)\^f{y)f{w). 

□ 

Now, setting / := Xi for a fixed index i G m} we note that there exists 

a vector subspace 14 in Hk{M) such that dim 14 = Hk{M) — such that 

Tf = f and Tj = p. We can take 14 to be the space spanned by the restrictions 
to M of all polynomials of total degree at most k — 1, i.e. Hk-i- The dimension 
of Nk = dim(iffc) is the Hilbert polynomial of degree for k > k^. Thus Nk = 
dfc" -I- where d is the degree of the variety M and n is the dimension. In 
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particular, denoting by Wk the orthogonal complement of 14 in Hk{M)r\L?{M, ^k) 
then dim(ll^fc) = Setting := Tj — Tf 2 gives 4^ = 0 on 14 and hence 

TrT| - TrT/2 = 0 + 

using that {TfPk,Pk) ^ sup |/|m (p/c,P/c)i2(M and dim 114 = 0{k^-^). 

□ 

3. Sampling and interpolation of real orthogonal polynomials 

3.1. Sampling polynomials in a real variety. 

Proof of Theorem [7J We equip M with the distance function d induced by the Eu¬ 
clidean distance in K™, i.e. d{x,y) := \x — y\. We recall that the corresponding 
Wasserstein L^-distance on the space 'P{M) of all probability measures on M is 
defined as 


W{p,(7)=mi d{x,y)dp{x,y), 

P JJmxm 


iMxM 

where the infimum is taken among all probability measures such that the first 
marginal of p is /i and the second a. The Wasserstein distance metrizes the weak-* 
convergence. 

We rely on the fact that 

1 


Nk 


Bk{x) dpkix) iy{x), 


where the convergence is in the weak-* topology, see [3]. Thus the way to prove 
the inequality of the theorem is by proving that there are constants {caIagS^j 
0 < ca < 1 such that 

W{ak,Pk)^0 

where at = Saga*, ^ai^a, Pk = '^Bk{x). Instead of the standard Wasserstein 
distance we will use an alternative expression more convenient for our purpose that 
it is equivalent to it, see [18]: 


lT(p, cr) = inf // d{x,y)\dp{x,y)\, 

P JjMxM 


iMxM 

where the inf is taken among all complex measures p such that the first marginal 
of / is p and the second a. The difference is that p is not necessarily positive and 
even if we don’t require that a and u are probability measures it still metrizes the 
weak-* convergence. Any candidate p with the right marginals is called a transport 
plan. 

The transport plan pk that is convenient to estimate is: 

Pk{x,y) = ^ X! X 9x(.x)^^y== dpk{x). 

where Kk{X,x) is the reproducing kernel for A in the space Hk and {g\}xeSk is 

the canonical dual frame (see (TUI) to < I in Hk- The latter is a frame 

Iv^JAeA. 

because Afc is sampling. 

If we compute the marginals of pk we get on one hand: 




Nk 


Ac Aj 


gxjx) 

x/MX) 




and the other marginal is given by 


dl3k{y) = 


Nk 


(x)^^^= dpk{x) = -^Kk{x,x)dpk{x). 

x/Bk{X) Nk 
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In the last equality we have used that g\ is a dual frame of the normalized repro¬ 
ducing kernels. 

The fact that {^a} it is the canonical dual frame to the normalized reproducing 
kernels allows us to conclude that = (SA(a^), positive and smaller 

than one. This follows from the following well known fact: 

Claim. If {Xn}n is a frame in a Hilbert space H and {t/n}™ is the dual frame then 

{Xn ; Vn) ^ [Oj 1] ■ 

Proof. Let T be the frame operator, i.e: T{x) = Since {xn}n is a 

frame then T is bounded, self-adjoint and invertible. The definition of the dual 
frame is T(yn) = Xn- For any vector v € H we have 


= T(T ^v)=^{T ^V,Xn)Xn. 


In particular 


Xk 


— 'y '. {vk: Xn)^ 


and multiplying by pk at both sides we get 

{xk.Vk) = X] K2/fc,x„)p. 

n 

Therefore {xk,yk) > 0 and {xk,yk) > 0 unless Xk = 0. Moreover, 
{xk,yk) - |( 2 /fc,Xfc)p = ^ |(j/fc,x„)p > 0. 


Thus 

therefore {xk,yk) < 1 too. 


n^k 

{xk,yk){I - {xk,yk)) > 0 , 


□ 


Finally we need to estimate 

I = ff \x - y\\dpk\ < ^ f \>^- x\\Kk{\x)\ ^^!!^ dpkix). 
JJmxm lik ^ Jm \/Bk{X) 


Since H^Alb — 1 we can estimate 

/2 < J_ ^ ^ 


AG Afc 


> M 


BkiX) 


We would like to use the sampling inequality (HU, and obtain that 


(3.1) 


A^Afc 

-ff 

Xlk JJmxM 


2 \Kk{X,x)\^ 

-x\ — . ' dyk{x) < 

M tSk{A) 

\y - x\’^\Kk[y,x)\’^ dp,k[x)dpk{y)- 


This we cannot do immediately because the polynomial (x — y)Kk{y,x) (in the 
variable y) is of degree fc -|- 1 instead of k as required in (II.1|) . 

But we are assuming that {p, (jf) define spaces with reproducing kernels of mod¬ 
erate growth. Thus Bk+i — Bk in M. Therefore if {Afejfc is sampling for Hk{M) 
then {Afc+i}fc is sampling for Hk{M). Thus, it is harmless to assume that A^ is 
sampling both for Hk and for Hk+i and we have established (EU- Then, using 
Theorem [HI we obtain 

W{ak,Pk) = 0{1/Vk), 

as desired. □ 










SAMPLING REAL POLYNOMIALS 


21 


3.2. Interpolating polynomials in a real variety. 

Definition 5. A sequence of sets of points on M is said to be interpolating for 
Hk{M) if the family of normalized reproducing kernels 

■.= Kk{;X)/\\Kk{-,X)\\ 

for A G Afc, is a Riesz sequence in the Hilbert space Hk{M), i.e.: 

^ |cAp < II Y ca«:a|P <C Y V{cA}AeA, G £2 

A^Afc A^Afc A^Afc 

where we will assume that C can be taken independent of k. 


This property is equivalent to the Plancherel-Polya inequality: 

(3.2) WfeHkiM) 

tt. Kk[X,X) 


and the interpolation property, for any sequence of sets of values {c^^^jAeA^ there 
are functions fk G Hk such that /fc(A(^^) = with 


(3.3) 




|ca| 


A^Afc 


KkiX,X)' 


and again the constant C should not depend on k. 

The property that the collection {KAjAGAt is a frame in H}~{M) is a quantitative 
version of the fact that the normalized reproducing kernels span the whole space 
and the property that they are a Riesz sequence quantifies the fact that they are 
linearly independent. 


Proof of Theorem [5l Let Fk C Hk be the subspace spanned by 
Kx{x) = Kk{X,x)/y^Kk{X,X) VAGAfc. 

Denote by g\ the dual (biorthogonal) basis to k\ in Fk- We have clearly that 

• We can span any function in Fk in terms of ka, thus: 

Y i^\{x)g\{x) = ICk{x,x), 

A^Afc 

where ICk{x,y) is the reproducing kernel of the subspace Fk- 

• The norm of g\ is uniformly bounded since k\ was a uniform Riesz sequence. 

• 5 a(A) = x/Kk{X, A). This is due to the biorthogonality and the reproducing 
property. 

We are going to prove that the measure (Jk = ^ SagA^ and the measure 
Pk = ■^ICk{x,x)dp{x) are very close to each other: W{ak,f3k) 0. In this case 
then since ICk{x,x) < Kk{x,x) and -^Bkix) dp, —>■ dv^ where v is the normalized 
equilibrium measure on M, then limsupj, Uk < i'- 

In order to prove that Wk(<JkT Pk) —>■ 0 we use the transport plan: 

Pk{x,y) = ^ Y gx{x)Kx{x) dp{x) 

AGAfc 

It has the right marginals, Cfc and (dk and we can estimate the integral in the same 
way as in the proof of Theorem [T] 

W{(7k,l3k)< [[ \x - y\\dpk\ = 0{1/Vk). 

JJmxm 
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The only point that merits a clarification is that we need an inequality similar to 
(13.11) . i.e: 


1 


E/ 


1 


\y - x\‘^\Kkiy, x)\‘^ dfi{x)dfi{y). 


MxM 


This time this is true because is a uniformly separated sequence by Proposi¬ 
tion m and therefore, it is a Plancherel-Polya sequence, see Proposition [T31 □ 


3.3. Sampling in convex domains. We proceed with the proof of Theorem U] 
The only part that we need to proof are the estimates for the reproducing kernel 
(11.51) . If these are proved, then it follows that the measure has the Bernstein- 
Markov property (EH) and the kernel is of moderate growth EH, thus we can 
apply Theorem [T] 


Proof. We start by the case when 17 = B is the unit ball. We denote by Bq (x) the 
Bergman function which is the reproducing kernel Kk{x,x) evaluated at the diag¬ 
onal of the space of polynomials of total degree k endowed with the norm with 
respect to the standard volume form. To get a lower bound for Bq we consider the 
cube Q such that the ball is inside it and tangent to its faces. Clearly by the compar¬ 
ison principle of the Bergman functions Bb > Bq and Bq(x) ~ Bi(xi) ■ ■ ■ Bi(xn) 
where Bj is the one dimensional Bergman kernel associated to the interval. This is 
known to be, see ED p. 108]: 

k 

\/d{x) ’ 

This implies that for points x in the interval that joins the origin with the center 
of one of the faces of the cube Q we have 




Bq{x) ~ min 



Thus we have the lower bound for Br that we wanted. To get the upper bound we 
will work in dimension n = 2 for simplicity but a similar argument works in any 
dimension. Observe that the space of polynomials of degree smaller or equal than 
k is spanned by the functions coS'^(t),/f>’sin'’(t)}j=o,...,fe in polar coordinates 
in the interval [0,1] x [0, 27r] with the measure pdp in the first interval and dt in 
the second. Consider now the space of functions in the product interval such 
that it is spanned by {p^ cos™(a:), sin™(x)}j=o,...,fc m=o,...,fc- The space Bfe is 

bigger than the space of polynomials thus the Bergman function at the diagonal 
> Bs^kix). But Bfj^ is easier to analyze because it is a product space of 
two one-dimensional spaces: The space of one dimensional polynomials of degree 
smaller k with the norm pdp in the interval [0,1] and the space of trigonometric 
polynomials {sin^ (x), cos-^ (a;)}j=o,...fc with the measure dx in [0,27r]. The Bergman 
function of Hk is the product of the one-dimensional Bergman functions. The 
Bergman function corresponding to the trigonometric polynomials is constant by 
invariance under rotations and by dimensionality it must be 2fc -|- 1. The space of 
polynomials in p are a space of Jacobi polynomials and its Bergman function has 
been estimated, see EH P- 108]: 


k 

Vd{x) 


Bj{x) ~ min 


Vx > 1/2. 
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Thus finally when n = 2 we get 

Bb < min ( —^=,kA . 

\ Vd{x) J 

Similarly in higher dimension we get 

C un \ 

y'dix) ) 

Now for an arbitrary convex domain there is an r > 0 (small) and an i? > 0 (big) 
that depend only on the domain such that for any point x in the boundary of the 
domain, there is a ball B{y, r) inside the domain tangent at x and with center y in 
the normal direction to the boundary of the domain at x and a cube Q{R) tangent 
to the domain at x in the middle of a face of the cube and such that the domain is 
contained in the cube. Again by the comparison principle of the Bergman function 
we get 


min 



^ Bq{r){x) < Ba{x) < BB(y,r){x) < min 




3.4. Existence of interpolating and sampling sequences in the one-dimensional 
setting. We conclude the paper by recalling some classical facts which are special 
for the one dimensional setting. 

Let /r be a finite measure on M with compact set K and assume that /i has the 
Bernstein-Markov property with respect to K. By the classical Christoffel-Darboux 
formula there exists constants ak+i such that 


Kk{x,y) = ttk+i 


qk+i{x)qk{y) - qk{x)qk+i{y) 

x-y 


where qk+i is the k th orthogonal polynomial (with respect to y). Let := {xj^^} 
be the set of fc + 1 zeros of qk (which by classical results are indeed all distinct 
and contained in the support K of /i). Then Kk{x[^\xj^^) = 0 ii i ^ j, as follows 
immediately from the Christoffel-Darboux formula. Hence, normalizing Kk (•, Xj ^'^) 
yields an orthonormal base in Hk {M, y) and as a consequence the following “sam¬ 
pling equality” holds for any pk G Hk{M) : 


f \pkfdy = '^ 

Jm 


Bk{x) 


Pk{xf'^) 


and in particular the sequence A^ is both sampling and interpolation. Finally, 
recall that by classical results the normalized Dirac measure 5k on the zeros Kk has 
the same weak limit points as Bk/{k -|- l)/i. In particular, if y has the Bernstein- 
Markov property, then 5k —>■ yeq, which is thus consistent with Theorem [1] and 
Theorem O 
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